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Introduction 


Since  1941,  the  State  of  Montana  has  used  post-season  surveys  of  hunting  license  and 
permit  (LPT)  holders  to  estimate  wildlife  harvest  (Thompson  1947).  Results  of  these  surveys  are 
used  to  develop  hunting  season  regulations,  evaluate  and  develop  wildlife  management 
strategies,  for  wildlife  research,  for  hunt  planning  by  the  public,  as  well  as  other  applications. 

Over  the  years,  methods  used  to  conduct  the  surveys  have  gone  through  several  changes. 
Some  of  these  changes  have  resulted  from  major  program  reviews  (Gooch  1971,  Cada  1983, 
Cada  et  al.  1986),  and  some  have  resulted  from  the  increasing  complexity  of  hunting  regulations 
(ex.  Gooch  1971  vs.  Martin  2005).  Several  program  changes  have  also  resulted  from  changes  in 
staff  responsible  for  conducting  the  harvest  surveys. 

License  years  2004  and  2005  mark  the  beginning  of  a  major  change  in  staffing  as  well  as 
a  reorganization  of  the  harvest  survey  team  and  program.  For  the  previous  20+  years,  the  survey 
team  has  consisted  of  staff  from  the  Research  and  Technical  Services  Bureau  (RTS)  in  the 
Wildlife  Division.  This  staff  included  the  Bureau  supervisor,  a  programmer/  analyst,  and  a 
survey  coordinator  to  oversee  a  staff  of  approximately  60-80  survey  interviewers  that  collected 
the  data.  The  team  has  often  had  additional  database  programming  support,  and  in  recent  years 
this  support  has  come  from  non-agency  contractors.  Beginning  in  license  year  2004,  the  RTS 
programmer/  analyst  position  was  vacated  and  then  eliminated.  Duties  related  to  data  entry  and 
data  management,  previously  a  responsibility  of  this  position,  have  been  transferred  to  the  new 
Application  Development  Bureau  in  the  Information  Services  Division.  This  permits  increased 
focus  on  these  activities  and  will  lead  to  substantial  program  development  in  coming  years. 
Subsequently,  RTS  created  a  biometrician  position,  which  will  have  significant  time  devoted  to 
the  harvest  surveys  program.  This  position  will  focus  on  procedures  for  sampling  LPT  holders 
and  estimating  hunting  and  harvest-related  quantities,  which  will  add  increased  focus  on  these 
areas.  This  will  lead  to  program  development  related  to  statistical  issues  in  coming  years. 
Additionally,  lead  survey  interviewer  positions  were  created  to  add  more  structure  to  the  team  of 
interviewers  and  to  share  some  of  the  supervisory  responsibilities  held  by  the  survey  coordinator. 
Finally,  a  working  group  was  formed  to  oversee  and  influence  the  direction  of  the  harvest 
surveys  program.  This  working  group  consists  of  approximately  15  biologists,  researchers,  and 
other  FWP  personnel  that  make  use  of  the  harvest  surveys  estimates  and  represent  various 
regions  and  divisions  of  the  agency.  The  scope  of  the  group  is  broad,  but  includes  program 
directions  as  well  as  estimation  methods. 

This  reorganization  has  resulted  in  several  changes  that  affect  estimates  of  hunting-  and 
harvest-  related  quantities.  As  this  is  the  first  year  these  changes  were  implemented,  methods  are 
described  with  some  detail  in  this  report.  Despite  these  changes,  tables  of  estimates  generated 
from  the  survey  data  look  similar  to  those  received  by  FWP  staff  in  previous  years.  Those  tables 
are  presented  in  the  results  section. 
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Methods 


Parameters  estimated 

In  an  effort  to  define  the  parameters  of  the  LPT  holder  population  that  need  to  be 
estimated  during  the  annual  hunting  and  harvest  survey,  an  effort  was  made  to  explicitly  identify 
what  quantities  related  to  antelope  hunting  and  harvest  are  necessary  for  annual  use  by  FWP  staff 
and  the  public. 

Hunting-related  parameters 

Hunter  success.  Based  on  recommendations  of  the  working  group,  it  was  determined  that  this 
quantity  is  useful  primarily  at  the  statewide  level  for  public  interest.  It  is  not  usually 
necessary  for  hunting  season  regulation  changes  at  the  hunting  district  or  regional  levels. 
It  is  defined  as  the  ratio  of  the  number  of  LPT  holders  that  harvested  at  least  one  antelope 
during  the  hunting  season  to  the  number  of  LPT  holders  that  actively  participated  in 
hunting.  Estimates  for  both  residents  and  nonresidents  are  required. 

Number  of  hunters.  This  quantity  is  estimated  at  the  hunting  district,  regional,  and  statewide 
levels.  It  is  often  useful  for  hunting  season  regulation  changes,  for  example  changes  in 
total  harvest  may  be  due  to  fewer  or  more  hunters  rather  than  changes  in  animal 
population  sizes.  Further,  it  is  required  to  estimate  the  economic  impact  of  hunting  in 
different  areas,  as  these  calculations  rely  on  the  number  of  hunters  and  hunter  days.  This 
quantity  is  defined  as  the  number  of  LPT  holders  that  spent  at  least  a  portion  of  one  day 
hunting  in  the  area  of  interest.  Note  that  because  each  LPT  holder  can  hunt  in  more  than 
one  hunting  district,  summing  the  total  number  of  hunters  in  each  hunting  district  will 
sum  to  more  than  the  actual  number  of  hunters  in  a  given  region  or  statewide.  If  a  LPT 
holder  hunts  in  more  than  one  hunting  district  in  a  region,  that  LPT  holder  will  count  for 
only  one  hunter  at  the  regional  level.  Estimates  of  hunters  at  the  statewide  level  are  based 
on  asking  LPT  holders  whether  or  not  they  hunted  at  all  during  the  license  year.  Further, 
because  LPT  holders  commonly  possess  and  simultaneously  hunt  with  more  than  one 
valid  LPT  at  the  same  time  in  a  given  hunting  district,  it  is  not  possible  to  make  estimates 
of  the  number  of  hunters  with  a  given  LPT  type.  An  LPT  holder  that  hunts  with  more 
than  one  valid  LPT  at  the  same  time  is  only  one  hunter;  they  are  not  counted  as  a  one 
hunter  for  each  LPT  they  possess  at  the  time  they  are  hunting.  Because  hunters  cannot  be 
both  residents  and  nonresidents,  however,  it  is  possible  to  estimate  the  number  of  resident 
and  non-resident  hunters  at  the  hunting  district,  regional,  and  statewide  levels. 

Hunter  days.  This  quantity  is  also  estimated  at  the  hunting  district,  regional,  and  statewide  level 
for  the  same  reasons  that  the  number  of  hunters  is  estimated  at  these  levels.  A  hunter  day 
is  defined  as  an  LPT  holder  spending  at  least  a  portion  of  a  day  hunting  in  a  particular 
hunting  district.  Therefore,  this  quantity  is  defined  as  the  number  of  days  that  LPT 
holders  spent  hunting  in  a  particular  hunting  district,  in  hunting  districts  that  fall  within  a 
given  region,  or  in  all  hunting  districts  in  the  state.  Note  that  it  is  possible  and  relatively 
common  to  hunt  in  more  than  one  hunting  district  on  one  day,  and  FWP  has  traditionally 
questioned  hunters  about  the  number  of  days  they  spent  hunting  in  each  hunting  district 
during  the  license  year  rather  than  each  region  or  statewide.  For  this  reason,  estimates  of 
the  total  number  of  hunter  days  in  a  region  or  statewide  will  be  more  than  the  number  of 
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days  that  LPT  holders  actually  spent  hunting  in  the  region  or  in  the  state.  This  quantity 
might  be  more  aptly  termed  hunter-district-days.  Also,  for  the  same  reasons  described 
above  for  estimates  of  the  number  of  hunters,  it  is  not  possible  to  estimate  the  number  of 
hunter  days  for  each  LPT  that  FWP  issues.  Estimates  of  hunter  days  are  made  for 
residents  and  nonresidents  at  the  hunting  district,  regional,  and  statewide  levels. 

Days  per  hunter  ratio.  This  quantity  is  estimated  at  the  hunting  district,  regional,  and  statewide 
levels.  It  is  used  to  monitor  changes  in  hunter  effort  as  an  index  of  the  density  and 
distribution  of  animals  with  a  given  area.  The  quantity  is  defined  as  the  ratio  of  the 
number  of  days  hunted  to  the  number  of  hunters  for  a  given  area,  and  thus  requires  no 
additional  questions  beyond  those  already  asked  of  LPT  holders  during  the  survey 
interview.  Because  the  number  of  hunter  days  and  the  number  of  hunters  are  estimated 
for  residents  and  nonresidents,  estimates  of  the  days  per  hunter  ratio  are  also  estimated 
for  residents  and  nonresidents  only  (i.e.  not  for  each  LPT). 

Harvest-related  parameters 

Total  harvest.  This  quantity  is  fundamental  to  wildlife  management  decisions  at  the  hunting 
district,  regional,  and  statewide  levels.  It  is  defined  as  the  total  number  of  animals 
harvested  by  LPT  holders  during  the  license  year. 

Management  success.  This  quantity  is  estimated  at  the  hunting  district,  regional,  and  statewide 
levels.  It  is  useful  for  making  changes  to  quotas  of  licenses  and  permits  at  these  levels.  It 
is  defined  as  the  total  harvest  estimate  divided  by  the  total  number  issued  of  a  given  LPT. 

Harvest  by  age-sex  classes.  This  quantity  is  also  fundamental  to  wildlife  management  in 

Montana.  Estimates  of  the  total  number  of  bucks,  does,  and  fawns  are  needed  at  the 
hunting  district,  regional,  and  statewide  levels  on  an  annual  basis  to  determine  if 
adjustments  to  hunting  season  regulations  or  quotas  need  to  be  made  in  order  to  meet 
management  goals. 

Harvest  by  time  period.  This  quantity  is  sometimes  useful  for  creating  hunting  regulations  to 

ensure  desired  distribution  of  harvest  across  time  during  the  hunting  season  at  the  hunting 
district,  regional,  or  statewide  levels.  In  previous  years,  this  quantity  was  defined  as  the 
total  harvest  in  each  of  seven  time  periods,  including  prior  to  the  rifle  season,  during  each 
week  of  the  four-week  rifle  season,  after  the  rifle  season,  and  an  unknown  category  for 
reported  kills  with  missing  data  for  time  period  of  kill.  For  this  report,  the  categories  were 
lumped  into  four  classes,  including  prior  to  the  rifle  season,  during  the  rifle  season,  after 
the  rifle  season,  and  an  unknown  category  for  reported  kills  with  missing  data  for  time 
period  of  kill.  This  will  give  us  a  greater  ability  to  have  precise  estimates  of  harvest  at  the 
hunting  district  level  for  relatively  small  estimates  of  harvest  in  each  of  these  categories. 

Sampling 

Strata  formation  and  sampling  rates 

Because  harvest-related  statistics  are  desired  and  can  be  estimated  by  residency  for  each 

LPT,  residency  and  LPT  were  used  to  define  sampling  strata.  Additionally,  the  FWP  Automated 

Licensing  System  (ALS)  stores  all  LPT  and  residency  information  as  subcategories  of  a  variable 

called  Item  Type,  which  indicates  how  certain  opportunities  were  sold  to  the  public.  Because 
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ALS  is  the  source  of  the  LPT  holder  population,  and  Item  Type  is  difficult  to  remove  from 
descriptions  of  opportunities  for  different  combinations  LPT  and  residency  within  the  ALS 
system,  Item  Type  was  also  included  in  the  formation  of  sampling  strata.  Thus,  sampling  strata 
were  defined  as  distinct  combinations  of  residency,  LPT  type,  and  Item  Type.  This  resulted  in 
formation  of  233  sampling  strata. 

Sampling  rates  for  these  strata  were  determined  in  large  part  using  methods  from 
previous  years.  Although  rationale  was  not  given  in  previous  years  for  sampling  rates  of 
particular  LPT’s,  the  harvest  survey  team  decided  not  to  make  large  changes  to  the  program  prior 
to  the  hiring  of  a  biometrician.  However,  the  team  did  seek  consultation  from  a  Montana  State 
University  graduate  student  in  Statistics  to  conduct  the  sampling.  The  contract  statistician  noted 
that  in  order  to  use  the  normal  approximation  to  construct  confidence  intervals  around  point 
estimates  of  hunting  and  harvest  statistics,  the  following  assumptions  are  made: 

(1)  np  >  10 

(2)  n{\- p)  >  10 

(3)  20  n<N , 

where  n  =  stratum  sample  size, 

N  =  stratum  population  size,  and 

p  =  the  estimated  proportion  of  LPT  holders  with  a  particular  trait  (for  example 
hunters  that  harvested  an  animal  within  a  particular  hunting  district). 

Based  on  these  assumptions,  as  well  as  methods  used  in  previous  years,  the  following 
rules  were  used  to  determine  sample  sizes  for  each  stratum: 

(1)  For  strata  with  up  to  375  LPT  holders,  all  LPT  holders  were  sampled. 

(2)  For  strata  with  1000  or  more  LPT  holders,  55%  of  the  LPT  holders  were 
sampled. 

(3)  For  strata  of  sizes  between  375  and  1000  LPT  holders,  the  percent  of  LPT 
holders  sampled  was  determined  based  on  a  linear  decrease  between  100% 
and  55%  of  LPT  holders  sampled. 

The  accuracy  and  precision  goals  for  these  sampling  rates  were  not  articulated. 

Sampling  design 

A  random  sample  of  size  specified  by  sampling  rate  rules  was  drawn  from  each  sampling 
stratum.  However,  owners  of  one  LPT  can  also  own  other  LPT’s,  meaning  that  the  sampling 
frames  for  the  various  strata  were  not  mutually  exclusive.  For  example,  many  doe-fawn  licenses 
require  an  either  sex  license  as  a  prerequisite.  Survey  protocol  requires  that  interviewers  ask  LPT 
holders  about  all  antelope  hunting  and  harvest  activities  regardless  of  which  LPT  they  were 
sampled  for.  This  implies  that  when  a  LPT  holder  was  drawn  into  the  sample  for  a  particular 
LPT,  they  were  also  included  in  the  sample  for  every  other  LPT  that  they  owned.  If  each  LPT 
owned  by  a  LPT  holder  is  considered  part  of  a  network  that  is  defined  by  the  ID  number  of  a 
particular  LPT  holder,  then  this  sampling  procedure  can  be  termed  “stratified  network  sampling.” 
Because  of  this  sampling  procedure,  it  is  difficult  to  sample  exactly  at  the  pre-determined 


5 


sampling  rates.  The  size  of  the  sample  for  each  stratum  depends  on  the  LPT’s  owned  by  the  LPT 
holders  in  that  stratum  and  the  sampling  rates  for  other  strata  to  which  those  LPT  holders  may 
also  belong.  Thus,  this  procedure  usually  results  in  over-sampling  of  strata. 

Data  collection 

A  team  of  approximately  65  telephone  interviewers  that  were  supervised  by  the  survey 
coordinator  collected  the  hunting  and  harvest  data.  Interviewers  were  hired  and  trained  during 
November  and  December  of  2005.  The  majority  of  interviewers  were  hired  in  Bozeman,  and  a 
survey  calling  center  was  developed  in  the  Bozeman  FWP  Regional  office.  Five  lead 
interviewers  were  hired  in  Bozeman  to  structure  the  interviewing  staff  at  the  call  center. 
Additionally,  approximately  a  third  of  the  interviewers  were  hired  in  other  areas  of  the  state  to 
call  resident  LPT  holders  residing  in  their  local  calling  areas. 

Calls  were  made  to  the  sample  of  antelope  LPT  holders  during  January  through  March  of 
2006.  Calls  were  made  primarily  from  cell  phones  to  avoid  long-distance  charges,  though  calls  in 
local  calling  areas  were  made  from  landline  phones  to  avoid  cell  phone  minute  charges.  Sampled 
LPT  holders  were  called  up  to  three  times  in  attempt  to  conduct  the  interview.  If  no  contact  was 
made  at  this  point,  uncontacted  LPT  holders  were  treated  as  non-respondents  to  the  survey. 

When  a  contact  was  established  with  an  LPT  holder,  interviewers  asked  license  holders  a  set  of 
pre -defined  questions  about  antelope  hunting  and  harvest  activities  during  the  2005  license  year 
(Appendix  1).  Answers  to  these  questions  were  recorded  on  paper  datasheets.  During  and 
following  the  calling  period,  lead  interviewers  and  the  survey  coordinator  edited  completed 
datasheets  in  an  effort  to  monitor  the  quality  of  data  and  to  correct  obvious  mistakes. 

After  the  data  were  collected  from  LPT  holders,  datasheets  were  bundled  and  delivered  to 
a  data  entry  contractor.  The  contractor  entered  the  data  into  a  centralized  Oracle  database  using 
data  entry  forms  developed  by  FWP  programmer  analysts  using  Oracle  Forms. 

Data  quality  control  and  management 

Quality  control 

Once  data  were  entered  to  the  database,  additional  quality  control  was  performed  using 
SQL  in  the  Oracle  database.  This  included  eliminating  duplicate  records  and  enforcing  data 
integrity  constraints  (e.g.  only  allowing  valid  hunting  districts,  etc.).  In  many  cases  this  involved 
the  creation  of  unknown  categories  for  analysis  and  estimation,  for  example  an  unknown  age-sex 
class,  unknown  hunting  districts  within  regions  and  statewide,  and  others.  While  these  categories 
may  not  be  of  direct  interest,  these  data  are  useful  for  estimating  hunting  and  harvest  parameters 
at  higher  levels,  for  example  total  harvest  including  the  unknown  age-sex  class  and  regional  and 
statewide  hunting  and  harvest  statistics  in  the  case  of  unknown  hunting  districts. 

Analysis  table  creation 


Because  hunting  effort  statistics  can  be  estimated  by  residency  status,  while  harvest 
statistics  can  be  estimated  by  residency  and  LPT,  two  tables  were  needed  for  analyses.  A  table  to 
estimate  hunting  effort  statistics  was  created  containing  a  record  for  every  resident  and  non¬ 
resident  license  holder  that  responded  to  the  antelope  survey.  Each  record  captures  information 
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on  whether  a  license  holder  hunted,  the  districts  in  which  they  hunted,  the  number  of  days  they 
hunted  in  those  districts,  and  whether  or  not  they  harvested  an  antelope. 

A  table  to  estimate  harvest  statistics  was  created  containing  a  record  for  every  LPT 
owned  by  every  license  holder  that  was  selected  for  the  antelope  survey.  Each  record  includes 
information  on  the  Item  Type  of  the  LPT,  which  hunter  owned  the  LPT,  the  residency  of  that 
hunter,  the  total  number  of  that  type  of  LPT  issued,  whether  the  hunter  that  owned  that  LPT 
responded  to  the  survey  or  not,  as  well  as  the  following  information  for  hunters  that  responded  to 
the  survey:  whether  an  antelope  was  harvested  with  each  LPT  owned  by  the  hunter,  as  well  as 
details  about  antelope  that  were  harvested  including  the  hunting  district  where  the  antelope  was 
harvested,  the  age-sex  class,  and  the  time  period  when  the  animal  was  harvested. 

A  major  issue  in  the  creation  of  the  table  for  estimation  of  harvest  statistics  is  that  we  did 
not  consistently  ask  LPT  holders  which  LPT  was  used  to  harvest  antelope  when  hunters  were 
successful.  Therefore  assigning  harvested  antelope  to  LPT’s  owned  by  a  successful  hunter 
required  additional  programming  work.  The  basic  algorithm  used  to  assign  harvested  antelope  to 
particular  LPT’s  relied  on  the  regulations  database.  Considering  the  legal  definitions  for  each 
LPT  owned  by  responding  LPT  holders,  such  as  allowable  age-sex  classes,  each  reported 
antelope  that  was  harvested  was  assigned  to  an  LPT  owned  by  the  LPT  holder. 

It  was  not  uncommon  for  a  harvested  animal  reported  by  a  hunter  to  fit  within  the  legal 
definitions  of  more  than  one  LPT,  for  example  some  hunting  districts  allow  hunters  to  harvest  a 
doe  antelope  with  an  either  sex  license  or  a  doe/fawn  license.  When  hunters  harvested  an  animal 
that  fit  within  the  legal  definitions  of  more  than  one  LPT  owned  by  the  hunter,  we  assumed  that 
hunters  harvested  antelope  with  LPT’s  in  order  of  most  restrictive  to  least  restrictive  hunting 
opportunities.  This  implies  that  hunters  made  use  of  doe/fawn  licenses  before  either  sex  licenses, 
that  they  used  LPT’s  valid  in  a  single  district  prior  to  LPT’s  valid  in  multiple  districts,  and  that 
they  used  LPT’s  limited  to  archery  equipment  before  LPT’s  with  no  limitations  on  method  of 
take.  The  primary  reason  for  making  this  assumption  is  that  if  hunters  own  more  than  one 
license,  using  LPT’s  to  harvest  animals  in  this  order  allows  hunters  maximum  flexibility  in 
remaining  hunting  opportunities  following  use  of  a  particular  LPT.  Bias  introduced  into 
estimates  of  harvest  by  making  this  assumption  would  be  an  over-estimate  of  harvest  with  more 
restrictive  LPT’s  and  an  underestimate  of  harvest  with  less  restrictive  LPT’s.  This  method  of 
assigning  kills  to  LPT’s  is  similar  to  that  used  in  previous  years. 

Another  uncertainty  that  arose  during  the  process  of  assigning  harvested  antelope  to 
LPT’s  owned  by  successful  hunters  was  the  case  that  characteristics  of  harvested  animals  did  not 
fit  within  the  legal  definitions  of  any  LPT(s)  owned  by  the  hunter.  The  most  common  of  these 
situations  were  when  a  hunter  owned  an  either  sex  license  and  a  doe/fawn  license  and  reported 
harvesting  two  antelope  bucks,  and  when  a  hunter  reported  harvesting  an  antelope  in  districts  for 
which  they  did  not  own  a  valid  LPT.  In  these  situations  (and  other  situations),  after  assigning 
reported  kills  to  LPT’s  for  which  those  kills  were  clearly  legal,  we  assigned  the  remaining 
reported  kills  to  the  remaining  LPT’s  regardless  of  whether  it  was  legal  for  a  hunter  to  harvest 
the  animals  with  those  LPT’s.  For  example  if  a  person  reports  killing  2  antelope  bucks  and  owns 
an  either  sex  license  and  a  doe/fawn  license,  we  assigned  one  kill  to  each  license.  This  differs 
from  methods  used  in  previous  years,  where  one  of  the  antelope  buck  kills  would  be  assigned  to 
the  either  sex  license,  and  the  other  would  be  assigned  to  the  "unknown"  LPT  (while  the 
doe/fawn  license  would  remain  unused).  We  chose  not  to  create  a  class  of  LPT  named 
“unknown.” 
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The  final  issue  that  arose  while  assigning  reported  antelope  kills  to  LPT’s  owned  by  a 
successful  hunter  was  the  case  where  people  reported  harvesting  more  antelope  than  they  had 
LPT  records  for  in  the  ALS  system.  In  these  cases,  these  extra  reported  kills  were  assigned  to  a 
“none”  LPT,  signifying  that  they  were  unlicensed  antelope  kills  made  by  hunters  that  had  at  least 
one  antelope  license.  Estimates  of  total  unlicensed  antelope  harvest  were  then  accordingly  made 
for  residents  and  nonresidents.  In  previous  years,  this  extra  antelope  would  be  assigned  to  the 
"unknown"  LPT  category  described  above. 

Data  analysis 

Response  rate 

Non-response  can  bias  estimates  of  hunting  effort  and  harvest  parameters.  The  potential 
for  bias  increases  with  the  amount  of  non-response,  and  therefore  higher  response  rates  protect 
against  non-response  bias.  In  an  effort  to  be  clear  about  the  potential  for  non-response  bias  in 
estimates  of  hunting  effort  and  harvest  statistics,  response  rates  were  calculated.  Because  hunting 
effort  statistics  were  calculated  for  residents  and  nonresidents  only,  while  harvest-related 
statistics  were  calculated  for  each  stratum,  response  rates  were  calculated  at  these  levels  as  the 
number  of  respondents  divided  by  the  number  of  LPT  holders  sampled  at  each  of  these  levels. 

Sampling  design  issues 

Analyzing  the  data  collected  from  each  stratum  as  though  it  were  a  simple  random 
sample  could  lead  to  biased  estimates  due  to  the  stratified  network  sampling  design.  For 
example,  a  fraction  of  the  nonresident  900-00  license  holder  population  is  composed  of  900-30 
license  holders,  because  the  900-00  license  is  a  prerequisite  for  900-30  licenses.  Because  the 
nonresident  900-30  population  is  less  than  375  individuals,  it  is  sampled  at  100%.  Thus, 
nonresident  900-30  holders  have  a  100%  chance  of  being  included  in  the  nonresident  900-00 
sample  because  of  the  stratified  network  sampling  design.  By  contrast,  nonresidents  that  only 
have  a  900-00  license  have  a  <100%  chance  of  being  included  in  the  sample  of  the  nonresident 
900-00  holder  population.  If  a  relationship  exists  between  the  probability  of  an  LPT  holder 
harvesting  an  antelope  with  a  particular  LPT  and  whether  or  not  that  LPT  holder  also  owns 
another  LPT,  the  estimate  of  total  harvest  for  the  LPT  holder  population  in  question  will  be 
biased  if  analyzed  as  though  it  were  a  simple  random  sample.  This  bias  will  be  proportional  to 
the  covariance  between  the  harvest  rate  with  the  LPT  in  question  and  the  harvest  rate  with  other 
LPT’s  for  hunters  that  have  two  types  of  LPT’s.  Specifically, 

bias(HB)*^U ’-HA  (1) 

ha 


where  H A  =  the  estimated  total  harvest  with  LPT  “A", 

Hb  =  the  true  harvest  rate  for  LPT  “B”,  and 
H A  =  the  true  harvest  rate  for  the  LPT  “A”  ( sensu  Bethlehem  2002). 
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To  protect  against  this  source  of  bias,  estimators  for  hunting  effort  and  harvest  must  be 
adjusted  to  account  for  variable  probabilities  of  inclusion  into  the  sample  for  different  sample 
networks.  These  inclusion  probabilities  can  be  defined  as  the  probability  that  any  LPT  attached 
to  a  particular  network  was  selected  in  any  stratum: 


nL 


=i-n 


y 


h=\ 


(Mh\ 


j 


(2) 


where  nk  =  the  inclusion  probability  for  LPT’s  in  the  network  owned  by  LPT  holder  k, 
h  =  stratum  number  1  to  L, 

Mh  =  the  number  of  LPT’s  issued  in  stratum  h, 

m*hk  =  the  number  of  LPT’s  attached  to  the  kth  network  in  stratum  h,  and 
nh  =  the  sample  size  in  stratum  h  (Thompson  2002). 


Inclusion  probabilities  were  calculated  only  for  responding  LPT  holders,  and  thus  we 
implicitly  assume  that  non-respondents  were  missing  at  random  (MAR)  from  the  sample  (Lohr 
1999).  For  this  method  to  produce  unbiased  estimates,  there  must  be  little  or  no  relationship 
between  the  probability  that  an  LPT  holder  responded  to  the  survey  and  the  hunting  effort  and 
harvest  statistics  being  estimated.  For  example,  if  non-respondents  were  less  likely  to  harvest  an 
antelope  than  respondents,  then  estimates  of  antelope  total  harvest  will  be  biased  high. 

Inclusion  probabilities  were  applied  to  each  responding  LPT  holder  in  the  case  of  hunting 
effort  related  statistics  and  to  every  LPT  owned  by  an  LPT  holder  in  the  case  of  harvest  related 
statistics  (Thompson  2002). 

Estimation  of  hunting  effort  and  harvest  parameters 

To  estimate  hunting  effort  statistics,  a  sampling  weight  was  calculated  for  each 
responding  LPT  holder  as: 


wk=—  (3) 

nk  ■ 

where  wk  =  the  sampling  weight  for  the  LPT  holder  that  owns  network  k,  and 
nk  =  the  inclusion  probability  for  network  k. 

To  use  these  weights  to  make  estimates  of  population  parameters,  they  should  reflect  the 
characteristics  of  the  population  of  interest  (Horvitz  and  Thompson  1952,  Lohr  1999).  In  this 
case,  this  implies  that  the  weights  for  responding  resident  LPT  holders  should  sum  to  the  total 
number  of  resident  LPT  holders  and  the  weights  for  responding  nonresident  LPT  holders  should 
sum  to  the  total  number  of  nonresident  LPT  holders.  Weights  were  therefore  “grossed  up”  to  the 
populations  of  resident  and  nonresident  LPT  holders  using: 
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(4) 


Nh 


where  wsk  =  the  grossed  up  weight  for  the  LPT  holder  that  owns  network  k, 
wk  =  the  sampling  weight  for  the  LPT  holder  that  owns  network  k, 

Nh  =  the  total  number  of  LPT  holders  in  residency  h,  and 

nh  =  the  number  of  responding  LPT  holders  in  residency  h  (Raab  et  al.  2006). 

Because  of  this  weighting  scheme,  calculations  to  estimate  hunting  effort  statistics 
implicitly  assume  that  non-responding  LPT  holders  are  MAR  within  each  residency.  Using 
grossed  up  weights,  hunter  success  was  estimated  at  the  statewide  level  as  a  weighted  mean 
success  rate  for  the  subpopulation  of  LPT  holders  that  reported  hunting: 


(5) 


where  sh  -  the  estimated  hunter  success  rate  for  residency  h, 

nh  =  the  number  of  responding  LPT  holders  in  residency  h  that  reported  hunting, 
yk  -  an  indicator  variable  indicating  whether  LPT  holder  k  harvested  one  or  more 
antelope,  and 

wsk  =  the  grossed  up  weight  for  the  LPT  holder  that  owns  network  k. 

The  number  of  hunters  was  estimated  at  the  statewide,  regional,  and  hunting  district 
levels  using  the  Horvitz-Thompson  estimator: 


(6) 


where  Hh  =  the  estimated  number  of  hunters  of  residency  h  in  the  area  of  interest, 
nh  =  the  number  of  responding  LPT  holders  in  residency  h, 
yk  =  a  variable  indicating  whether  LPT  holder  k  hunted  in  the  area  of  interest,  and 
wsk  =  the  grossed  up  weight  for  the  LPT  holder  that  owns  network  k  (Horvitz  and 


Thompson  1952). 


To  estimate  the  number  of  hunter  days  and  the  days:hunter  ratio  at  the  statewide, 
regional,  and  hunting  district  levels,  we  first  had  to  deal  with  the  issue  of  item  non-response  for 
the  number  of  days  reported  hunting.  That  is,  7  hunters  reported  hunting  in  particular  district  in 
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license  year  2005  without  reporting  the  number  of  days  that  they  hunted  in  that  district.  In  order 
to  use  the  full  sample  to  estimate  the  number  of  hunter  days  and  the  days:hunter  ratio,  we 
decided  to  impute  data  for  the  number  of  days  hunted  by  those  7  hunters  that  did  not  report  the 
number  of  days  that  they  hunted.  In  order  to  eliminate  effects  on  the  estimated  total  number  of 
days  hunted  and  the  days:hunter  ratio,  we  could  set  the  value  of  the  missing  data  equal  to  the 
mean  for  residents  and  nonresidents.  This  “cell  mean  imputation”  would  effectively  not  use  the 
imputed  data  values  in  estimating  the  total  number  of  days  or  the  days:hunter  ratio.  However, 
this  would  lead  to  underestimation  of  our  uncertainty  in  those  estimates  because  these  7  hunters 
would  be  assigned  days  hunted  values  that  did  not  deviate  from  the  mean.  In  order  to  minimize 
effects  on  estimates  of  the  total  number  of  days  hunted  and  the  days:hunter  ratio,  as  well  as  the 
associated  estimates  of  uncertainty,  we  therefore  use  a  “cell  mean  model”  to  impute  these 
missing  values  (Lohr  1999).  We  assumed  that  these  hunters  had  days  hunted  values  that  fell 
along  a  normal  distribution  with  mean  equal  to  the  weighted  mean  for  residents  or  nonresidents, 
and  variance  equal  to  the  weighted  variance  for  residents  or  nonresidents: 

d'u,  ~N(Mk,el)  (7) 

where  d*hk  =  the  missing  value  for  days  hunted  for  the  LPT  of  residency  h  that  owns 
network  k , 

fxh  =  the  weighted  mean  number  of  days  hunted  for  residency  h,  and 

=  the  weighted  variance  of  the  mean  number  of  days  hunted  for  residency  h. 

A 

Once  these  7  values  were  imputed,  the  number  of  days  hunted  for  residency  h{Dh )  was 

estimated  using  the  Horvitz-Thompson  estimator  in  equation  (6)  above,  with  yk  set  equal  to  the 

number  of  days  hunted  for  LPT  holder  k  in  the  hunting  district  of  interest  (or  the  number  of  days 
hunted  in  hunting  districts  falling  within  the  region  of  interested  or  the  total  number  of  days 
hunted  in  hunting  districts  within  the  entire  state  at  the  higher  levels).  The  days:hunter  ratio  for 
each  area  of  interest  was  then  estimated  using  a  ratio  estimator: 


R, 


(8) 


A 

where  Rh=  the  estimated  days:hunter  ratio  in  the  area  of  interest  for  residency  h, 

A 

Dh  =  the  estimated  number  of  hunter  days  in  the  area  of  interest  for  residency  h, 
and 

A 

Hh  =  the  estimated  number  of  hunters  in  the  area  of  interest  for  residency  h. 


All  variances  for  hunter  effort  statistics  were  estimated  using  Taylor  Series  linearization 
(Woodruff  1971)  and  corrected  for  large  sampling  fractions  using  a  finite  population  correction 
(fpc)  factor  (Thompson  2002).  Combined  residency  hunter  effort  statistics  were  estimated  using 
equations  for  stratified  sampling,  treating  residency  as  the  stratifying  variable  (Thompson  2002). 
To  estimate  the  combined  residency  days:hunter  ratio  in  a  particular  area,  estimates  for  residents 
and  nonresidents  were  weighted  by  the  estimated  number  of  hunters  in  that  area  rather  than  the 


11 


total  number  of  resident  and  nonresident  LPT  holders.  This  was  done  in  order  to  avoid  down¬ 
weighting  estimates  of  the  ratio  in  the  cases  where  only  residents  or  nonresidents  reported 
hunting  in  an  area.  Therefore,  combined  residency  estimates  of  the  days:hunter  ratio  are 
conditioned  on  the  LPT  holders  that  reported  hunting  in  a  particular  area,  and  can  be  interpreted 
as  the  expected  number  of  days  that  an  LPT  holder  hunted  in  a  particular  area,  given  that  they 
hunted  in  the  area. 

Confidence  intervals  around  hunter  effort  statistics  were  calculated  using  the  normal 
distribution  approximation: 


y  ±  z*  yl\ar(y)  (9) 

where  y  =  the  point  estimate  for  the  population  quantity  of  interest, 

the  90th  percentile  value  from  the  normal  distribution  (1.2816),  and 
var(y)  =  the  fpc-corrected  sample  variance  of  the  point  estimate  for  the 
population  quantity  of  interest. 

The  absolute  lower  boundary  of  confidence  intervals  was  set  to  the  number  of  LPT  holders  in  the 
sample  with  the  population  characteristic  of  interest,  as  this  is  the  minimum  number  of  LPT 
holders  known  to  have  the  characteristic.  Confidence  intervals  were  not  adjusted  for  multiple 
comparisons,  and  therefore  cannot  be  interpreted  as  simultaneous  confidence  intervals.  Thus  we 
expect  20%  of  the  intervals  around  point  estimates  not  to  contain  the  true  population  value. 

To  estimate  antelope  harvest  statistics,  the  sampling  weight  for  each  responding  LPT 
holder  was  applied  to  every  LPT  owned  by  that  LPT  holder  (Thompson  2002).  These  LPT 
sampling  weights  were  then  grossed  up  to  sum  to  the  number  of  LPT’s  issued  within  each 
sampling  stratum  using  formula  (4)  above.  This  weighting  scheme  assumes  that  LPT’s  owned  by 
non-responding  LPT  holders  were  MAR  within  each  sampling  stratum. 

Antelope  harvest  statistics  were  estimated  as  totals  for  hunting  districts,  regions,  and  the 
state  by  residency  and  LPT  or  residency  and  LPT  types  (ex.  either-sex,  doe/fawn,  etc.).  Total 
harvest,  as  well  as  harvest  by  age-sex  classes,  and  time  period  within  each  stratum  in  the  area  of 
interest  were  estimated  with  the  Horvitz-Thompson  estimator  in  equation  (6)  above.  Estimates 
were  generated  for  each  sampling  stratum,  and  then  combined  into  estimates  for  residency  and 
LPT,  by  LPT  across  residency,  and/or  by  type  of  LPT  using  estimators  for  stratified  sampling 
(Thompson  2002). 

To  estimate  harvest  made  without  LPT’s,  a  separate  table  was  created  for  reported 
harvest  linked  to  the  “none”  LPT  for  hunters  that  reported  harvesting  more  antelope  than  the 
number  of  LPT’s  they  owned.  A  record  was  then  entered  in  this  table  for  every  resident  and 
nonresident  that  responded  to  the  survey  and  did  not  report  harvesting  more  antelope  than  the 
number  of  LPT’s  they  owned.  Sampling  weights  for  each  hunter  in  the  table  were  then  grossed 
up  to  sum  to  the  number  of  LPT’s  issued  within  each  residency  using  formula  (4)  above.  This 
new  table  was  then  used  to  estimate  unlicensed  total  harvest,  harvest  by  age-sex  classes,  and 
harvest  by  time  period  for  residents  and  non-residents  that  purchased  at  least  one  license.  These 
estimates  were  then  combined  with  estimates  of  harvest  for  antelope  legally  harvested  with 
LPT’s  using  formulas  for  stratified  sampling  (Thompson  2002). 

Following  tradition  from  the  previous  decade,  the  only  harvest  statistic  for  which  an  80% 
confidence  interval  was  estimated  was  total  harvest.  The  variance  of  total  harvest  estimates  was 
estimated  using  Taylor  Series  linearization  (Woodruff  1971)  and  was  fpc  corrected  (Thompson 


12 


2002).  However,  the  responding  sample  size  for  some  strata  was  only  one  individual.  In  these 
cases,  the  stratum  contribution  to  the  variance  is  taken  to  be  the  average  of  all  the  strata  with  a 
responding  sample  size  of  more  than  one  individual.  This  adjustment  results  in  a  conservative 
estimate  of  the  total  variance  (Lumley  2004).  Estimates  of  population  parameters  and  associated 
variances  were  not  made  for  strata  in  which  there  were  no  respondents;  all  statistics  for  these 
strata  were  set  equal  to  zero. 

The  80%  confidence  interval  on  total  harvest  was  estimated  using  the  normal 
approximation  in  equation  (9)  above.  The  absolute  lower  boundary  of  the  total  harvest 
confidence  intervals  were  set  to  the  number  of  LPT  holders  that  reported  harvesting  a  antelope  in 
the  area  of  interest,  as  we  know  that  this  is  the  minimum  number  of  antelope  harvested  in  each 
area  of  interest. 

All  data  analyses  were  conducted  using  the  R  statistical  software  program  (R 
Development  Core  Team  2006).  Extensive  use  was  made  of  the  survey  analysis  package  written 
for  R  (Lumley  2004). 

Results 


Response  rates  were  generated  for  each  LPT  that  was  issued  for  antelope  hunting 
privileges,  as  well  as  by  residency  (Table  1).  Estimates  of  hunter  success  were  generated  only  at 
the  statewide  level  (Table  2).  Following  requests  from  the  harvest  survey  working  group,  the 
remaining  hunter  effort  statistics  were  included  in  the  same  table  containing  estimates  of  total 
harvest,  harvest  by  age-sex  classes,  and  harvest  by  time  period  (Table  3). 
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Table  1.  Response  rates  for  antelope  hunting  and  harvest  survey  strata,  license  year  2005. 


residency 

LPT 

number  issued 

number sampled 

responses 

response  rate 

N 

300-00 

26 

26 

16 

0.62 

N 

300-10 

6 

6 

4 

0.67 

N 

301-00 

7 

7 

5 

0.71 

N 

310-00 

3 

3 

3 

1 

N 

311-00 

1 

1 

0 

0 

N 

318-00 

11 

11 

8 

0.73 

N 

320-00 

9 

9 

7 

0.78 

N 

320-10 

9 

9 

6 

0.67 

N 

321-00 

9 

9 

8 

0.89 

N 

321-10 

2 

2 

1 

0.5 

N 

329-00 

16 

16 

11 

0.69 

N 

329-10 

3 

3 

1 

0.33 

N 

330-00 

36 

36 

24 

0.67 

N 

330-10 

6 

6 

4 

0.67 

N 

338-00 

38 

38 

24 

0.63 

N 

339-00 

10 

10 

2 

0.2 

N 

340-00 

3 

3 

2 

0.67 

N 

341-00 

7 

7 

3 

0.43 

N 

350-00 

1 

1 

1 

1 

N 

360-00 

46 

46 

30 

0.65 

N 

360-03 

5 

5 

4 

0.8 

N 

370-10 

2 

2 

2 

1 

N 

380-00 

1 

1 

0 

0 

N 

380-10 

2 

2 

1 

0.5 

N 

390-00 

9 

9 

8 

0.89 

N 

401-00 

14 

14 

10 

0.71 

N 

404-00 

7 

7 

6 

0.86 

N 

413-00 

10 

10 

5 

0.5 

N 

413-10 

2 

2 

2 

1 

N 

420-00 

40 

40 

28 

0.7 

N 

420-10 

90 

90 

57 

0.63 

N 

430-00 

30 

30 

20 

0.67 

N 

430-10 

25 

25 

19 

0.76 

N 

440-00 

1 

1 

0 

0 

N 

440-10 

1 

1 

1 

1 

N 

444-00 

2 

2 

1 

0.5 

N 

444-10 

5 

5 

3 

0.6 

N 

450-00 

30 

30 

24 

0.8 

N 

450-10 

41 

41 

26 

0.63 

N 

470-00 

21 

21 

14 

0.67 

N 

470-10 

2 

2 

2 

1 

N 

471-00 

10 

10 

7 

0.7 

N 

471-10 

21 

21 

15 

0.71 
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Table  1.  Response  rates  for  antelope  hunting  and  harvest  survey  strata,  license  year  2005. 


residency 

LPT 

numberjssued 

number  sampled 

responses 

response  rate 

N 

473-00 

4 

4 

2 

0.5 

N 

480-00 

113 

113 

74 

0.65 

N 

480-10 

199 

199 

125 

0.63 

N 

481-00 

118 

118 

83 

0.7 

N 

481-10 

148 

148 

107 

0.72 

N 

490-00 

34 

34 

28 

0.82 

N 

490-10 

60 

60 

33 

0.55 

N 

491-00 

17 

17 

11 

0.65 

N 

491-10 

22 

22 

11 

0.5 

N 

500-00 

47 

47 

34 

0.72 

N 

500-30 

38 

38 

25 

0.66 

N 

501-00 

69 

69 

45 

0.65 

N 

501-30 

25 

25 

16 

0.64 

N 

510-00 

40 

40 

31 

0.78 

N 

511-00 

11 

11 

9 

0.82 

N 

512-00 

7 

7 

4 

0.57 

N 

513-00 

98 

98 

62 

0.63 

N 

513-30 

52 

52 

27 

0.52 

N 

514-00 

4 

4 

2 

0.5 

N 

530-00 

150 

148 

104 

0.7 

N 

530-30 

75 

75 

51 

0.68 

N 

540-00 

28 

28 

21 

0.75 

N 

550-00 

8 

8 

4 

0.5 

N 

560-00 

70 

70 

42 

0.6 

N 

570-00 

60 

60 

37 

0.62 

N 

571-00 

40 

40 

32 

0.8 

N 

572-00 

24 

24 

18 

0.75 

N 

590-00 

67 

67 

42 

0.63 

N 

590-30 

37 

37 

21 

0.57 

N 

600-00 

49 

49 

27 

0.55 

N 

600-10 

17 

17 

13 

0.76 

N 

620-00 

43 

43 

25 

0.58 

N 

620-10 

10 

10 

8 

0.8 

N 

630-00 

13 

13 

8 

0.62 

N 

630-10 

1 

1 

1 

1 

N 

650-00 

40 

40 

28 

0.7 

N 

650-10 

10 

10 

7 

0.7 

N 

670-00 

34 

34 

24 

0.71 

N 
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3 

3 

1 

0.33 

N 
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98 

98 
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0.74 
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33 

0.65 

N 
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0.67 

N 
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1016 
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Table  1.  Response  rates  for  antelope  hunting  and  harvest  survey  strata,  license  year  2005. 


residency 

LPT 

numberjssued 

number  sampled 

responses 

response  rate 

N 

900-00 

811 

738 

453 

0.61 

N 

900-30 

89 

89 

56 

0.63 

N 

Archery Only 

900 

827 

509 

0.62 

N 

Doe  Fawn 

3118 

2527 

1643 

0.65 

N 

Either  Sex 

2975 

2331 

1569 

0.67 

N 

SUM 

6993 

5685 

3721 

0.65 

R 

215-00 

10 

10 

6 

0.6 

R 

215-10 

15 

15 

12 

0.8 

R 

300-00 

274 

274 

204 

0.74 

R 

300-10 

94 

94 

72 

0.77 

R 

301-00 

68 

68 

43 

0.63 

R 

310-00 

97 

97 

77 

0.79 

R 

310-10 

50 

50 

32 

0.64 

R 

311-00 

74 

74 

56 

0.76 

R 

318-00 

189 

189 

134 

0.71 

R 

319-00 

50 

50 

36 

0.72 

R 

320-00 

190 

190 

148 

0.78 

R 

320-10 

191 

191 

140 

0.73 

R 

321-00 

141 

141 

94 

0.67 

R 

321-10 

23 

23 

17 

0.74 

R 

329-00 

234 

234 

171 

0.73 

R 

329-10 

47 

47 

34 

0.72 

R 

330-00 

365 

365 

266 

0.73 

R 

330-10 

142 

142 

111 

0.78 

R 

338-00 

360 

360 

259 

0.72 

R 

339-00 

90 

90 

65 

0.72 

R 

340-00 

72 

72 

50 

0.69 

R 

341-00 

143 

143 

112 

0.78 

R 

350-00 

39 

39 

31 

0.79 

R 

350-10 

10 

10 

7 

0.7 

R 

360-00 

454 

379 

266 

0.7 

R 

360-03 

95 

95 

67 

0.71 

R 

370-00 

50 

50 

37 

0.74 

R 

370-10 

48 

48 

37 

0.77 

R 

371-00 

35 

35 

26 

0.74 

R 

371-10 

50 

50 

36 

0.72 

R 

380-00 

74 

74 

52 

0.7 

R 

380-01 

35 

35 

25 

0.71 

R 

380-10 

97 

97 

66 

0.68 

R 

381-00 

10 

10 

9 

0.9 

R 

381-10 

15 

15 

9 

0.6 

R 

381-11 

15 

15 

8 

0.53 

R 

390-00 

190 

190 

116 

0.61 
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Table  1.  Response  rates  for  antelope  hunting  and  harvest  survey  strata,  license  year  2005. 


residency 

LPT 

numberjssued 

number  sampled 

responses 

response  rate 

R 

401-00 

386 

375 

283 

0.75 

R 

404-00 

192 

192 

136 

0.71 

R 

413-00 

191 

191 

140 

0.73 

R 

413-10 

98 

98 

66 

0.67 

R 

420-00 

360 

360 

266 

0.74 

R 

420-10 

505 

474 

294 

0.62 

R 

430-00 

271 

271 

190 

0.7 

R 

430-10 

450 

450 

292 

0.65 

R 

440-00 

74 

74 

50 

0.68 

R 

440-10 

99 

99 

70 

0.71 

R 

441-00 

10 

10 

6 

0.6 

R 

444-00 

98 

98 

75 

0.77 

R 

444-10 

95 

95 

71 

0.75 

R 

450-00 

370 

370 

288 

0.78 

R 

450-10 

1158 

669 

496 

0.74 

R 

455-00 

2 

2 

2 

1 

R 

470-00 

229 

229 

152 

0.66 

R 

470-10 

48 

48 

37 

0.77 

R 

471-00 

190 

190 

147 

0.77 

R 

471-10 

329 

329 

221 

0.67 

R 

473-00 

146 

146 

97 

0.66 

R 

480-00 

1086 

476 

330 

0.69 

R 

480-10 

1197 

585 

399 

0.68 

R 

481-00 

1080 

485 

335 

0.69 

R 

481-10 

854 

470 

301 

0.64 

R 

490-00 

316 

316 

204 

0.65 

R 

490-10 

690 

546 

356 

0.65 

R 

491-00 

158 

158 

113 

0.72 

R 

491-10 

378 

378 

261 

0.69 

R 

500-00 

553 

382 

265 

0.69 

R 

500-30 

463 

449 

310 

0.69 

R 

501-00 

680 

392 

277 

0.71 

R 

501-30 

498 

461 

325 

0.7 

R 

510-00 

360 

360 

255 

0.71 

R 

511-00 

239 

239 

174 

0.73 

R 

512-00 

68 

68 

52 

0.76 

R 

513-00 

903 

401 

280 

0.7 

R 

513-30 

598 

522 

356 

0.68 

R 

514-00 

71 

71 

53 

0.75 

R 

530-00 

1353 

546 

391 

0.72 

R 

530-30 

916 

696 

505 

0.73 

R 

540-00 

269 

269 

204 

0.76 

R 

550-00 

292 

292 

202 

0.69 
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Table  1.  Response  rates  for  antelope  hunting  and  harvest  survey  strata,  license  year  2005. 


residency 

LPT 

number  issued 

number  sampled 

responses 

response rate 

R 

560-00 

931 

398 

288 

0.72 

R 

570-00 

540 

384 

272 

0.71 

R 

571-00 

359 

359 

273 

0.76 

R 

572-00 

225 

225 

184 

0.82 

R 

590-00 

632 

387 

277 

0.72 

R 

590-30 

414 

411 

294 

0.72 

R 

600-00 

855 

421 

301 

0.71 

R 

600-10 

208 

208 

148 

0.71 

R 

620-00 

456 

381 

268 

0.7 

R 

620-10 

90 

90 

72 

0.8 

R 

630-00 

137 

137 

113 

0.82 

R 

630-10 

9 

9 

9 

1 

R 

650-00 

360 

360 

268 

0.74 

R 

650-10 

38 

38 

28 

0.74 

R 

670-00 

316 

316 

219 

0.69 

R 

670-10 

12 

12 

8 

0.67 

R 

690-00 

903 

428 

322 

0.75 

R 

690-10 

472 

472 

337 

0.71 

R 

700-00 

11697 

4728 

3353 

0.71 

R 

700-20 

5846 

4300 

2994 

0.7 

R 

900-00 

3661 

1499 

1021 

0.68 

R 

900-30 

438 

409 

281 

0.69 

R 

Archery Only 

4099 

1908 

1302 

0.68 

R 

Doe  Fawn 

16342 

12786 

8890 

0.7 

R 

Either Sex 

30617 

18671 

13396 

0.72 

R 

SUM 

51058 

33365 

23588 

0.71 

SUM 

SUM 

58051 

39050 

27309 

0.7 
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Table  2.  Antelope  hunting  success  estimates  for  residents  and  nonresidents,  license  year  2005. 


district 

Res 

pet  success 

pet  succ  LCI80 

pet  succ  UCI80 

State 

N 

77.15 

76.37 

77.94 

State 

R 

65.63 

65.17 

66.1 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Table  3.  Antelope  hunting  effort,  total  harvest,  age-sex  class  harvest,  and  time  period  of  harvest  estimates,  license  year  2005. 
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Discussion 


License  year  2005  marks  the  beginning  of  several  changes  in  the  harvest  survey  program. 
One  of  these  changes  is  collapsing  the  many  antelope  time  period  of  harvest  classes  for  which 
estimates  were  made  in  previous  years  into  smaller  number  of  categories.  This  was  done  in  order 
to  make  more  reliable  estimates  of  harvest  for  these  parameters.  When  we  estimate  harvest  in 
these  categories  in  each  hunting  district  for  individual  LPT’s,  our  ability  to  meet  accuracy  and 
precision  goals  is  limited.  At  working  group  meetings  following  the  collection  of  these  data,  the 
desired  level  of  precision  was  stated  as  <  ±15%  of  point  estimates,  using  80%  confidence 
intervals.  Using  this  metric,  approximately  55%  of  the  total  harvest  point  estimates  contained  in 
the  table  of  estimates  above  do  not  meet  our  precision  goals.  Because  estimating  harvest  by  time 
period  categories  always  involves  estimating  smaller  quantities  than  when  estimating  total 
harvest,  our  ability  to  reliably  estimate  quantities  in  these  categories  will  be  even  more  limited. 

In  coming  years,  precision  goals  will  form  the  basis  of  sampling  procedures,  and  at  this  time  we 
will  formally  evaluate  our  ability  to  dependably  estimate  harvest  in  the  reduced  number  of  time 
period  categories. 

Other  changes  included  the  assignment  of  kills  to  LPT’s  which  could  not  legally  be  used 
to  make  the  described  kill  and  the  estimation  of  unlicensed  antelope  harvest  by  LPT  holders  that 
owned  at  least  one  license.  These  situations  have  apparently  existed  in  the  data  collected  for 
harvest  surveys  since  at  least  1995.  In  previous  years,  reported  kills  that  fell  into  these  categories 
were  assigned  to  a  single  “unknown”  LPT  category,  for  which  estimates  of  harvest  were  made. 
The  decision  to  change  this  practice  was  made  because  we  thought  it  would  be  useful  managers 
to  know  what  proportion  of  the  licenses  and  permits  they  issue  are  used  in  the  way  that  they 
intended,  as  well  as  how  many  were  not.  Accordingly,  we  thought  that  managers  could  make  use 
of  estimates  of  harvest  that  was  made  without  LPT’s. 

A  major  forthcoming  change  is  in  regards  to  the  assignment  of  harvested  antelope  to 
specific  LPT’s  owned  by  an  LPT  holder.  The  practice  of  not  asking  successful  hunters  which 
LPT  was  used  to  harvest  an  antelope  is  the  most  difficult  and  time  consuming  issue  to  deal  with 
prior  to  data  analyses.  It  also  undoubtedly  introduces  bias  that  would  be  difficult  to  quantify  and 
correct  for.  This  practice  will  be  discontinued  in  coming  years  in  favor  of  asking  LPT  holders 
about  harvest  success  and  characteristics  with  individual  LPT’s  that  they  own.  This  may  also 
eliminate  the  need  to  calculate  inclusion  probabilities  for  harvest-related  statistics,  because  we 
will  not  need  to  ask  LPT  holders  about  antelope  harvested  with  all  of  the  LPT’s  that  they  own. 
This  will  permit  harvest  estimation  similar  to  that  conducted  during  the  previous  decade,  when 
inclusion  probabilities  were  not  accounted  for.  However,  it  is  clear  that  variable  inclusion 
probabilities  existed  for  LPT  holders  during  this  time  frame,  so  the  potential  bias  described  in 
equation  (1)  above  is  implicit  in  hunting  and  harvest  estimates  dating  to  at  least  1995. 

Also,  because  confidence  intervals  for  harvest  statistics  are  only  estimated  around  total 
harvest,  biologists  that  use  the  various  harvest  estimates  for  season  setting  cannot  assess  the 
reliability  and  precision  of  most  estimated  harvest  parameters.  This  issue  will  be  addressed  in 
coming  years.  Estimates  that  do  not  meet  precision  goals  will  be  flagged  in  reports. 

Despite  the  changes  in  staff  and  methods,  presentation  of  hunting  and  harvest  estimates  is 
comparable  to  estimates  from  previous  years.  Many  program  changes  may  go  unnoticed  by  end- 
users  of  the  estimates  from  the  hunting  and  harvest  surveys.  This  will  in  fact  be  a  goal  of  the 
program  during  the  transition  period,  such  that  improvements  to  the  program  are  more  likely  to 
gain  acceptance  by  end-users. 
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Appendix  1:  2005  Antelope  hunting  and  harvest  survey  questionnaire. 
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'r  which  yee  Kilted  you:  ante  ops. _ ►  ' 

I  Sea  Anielap«  Hunting  Dislnti  V»p) 

3.  I  low  mar.y  amelcpo  d  6  you  ki=l  tr.  Montana  in  2005'? . 


.□ 


i .  ;  4<<xfex  .  •  Jtfyr  <  » :  v 


4.  h  which  antelope  ha  ring  cist'ict  did  you  !<::*  this  mtolapa?.... 


I 


AttTU.CiZ  » 
MOHNto! 


>1  • 

ft  i! 

y 

.h  :m<2  „ 


jl  5.  Antelope  kllieu  w»  a:  CJ^FAWN  Q* ADULT  DOE  Q3A Dll LT  BUCK 

n  fi.  An  I  dope  was  killed  during  wh  Oh  lima  per:ed7  u 

I  Q1  Before  Oct  9  OaQ«  906  QJOCt  17-23  Q40ct  24-  Ocf  30  O*0«t  31  •  No«  8  O®  After  Np«  S  \\ 
” :  •  >  i* 1  safe 'rnmm~i»tn  SEffTT??:  T? k  :*.y r  *z :  ••  ir: 

?n  whict*  antel0.dS  kuntrg  difctrct  file  you  kilt  th.s  ante-epe? . 

p  s  Asikiapc  Kltiod  woa  a:  OVawm  O2adultd0e  O3  ADULT  BUCK 

]<j  6.  Antelope  was  killed  during  which  ti.mepa'ioc? 


.^vs, 

i  i 

l 

rr« 

^  \ — >vl-«  j 

AVTElC^e  M.WTVWG 
OiS7riiC7 


J  O1  Before  OCt  9  O20ct  9  -  >6  O30ct  17-23  O40cf  24  -  Ott  30  O^Oci  31  ■  New  6  O®  Alter  NOV  6 

-;-7.  zs>-.'vr..-’xyr ■'■?•■  'mr 


i£ 


4,  Ip  Which  antelope  huntrn  distrsf  die  you  kill  Ihs  anlo:cpc? 


. rr 


.  5.  A'l^looo  killed  was  a: 


OVaWH  ■'J>?ADULTOOE  O3 ADULT  BUCK  AWrr. rWmYTO? 


(>.  Aateiopo  wa«  killed  during  which  time  period? 


fASjffsr 

period? 

p!:  O1  Bel  ore  Oct  9  O70cl  9-16  O50cf  17  -  23  C>4Oet  24  -  Oct  3D  O50ci  31  -  Nov  6  O^Altar  Nov  6 

ra:w,  w ■&. ^  ^  **.,*»■-* 
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Appendix  2 


Microsoft  Excel  file  containing  hunting  and  harvest  estimates,  available  electronically 
from  the  FWP  website  and  upon  request  from  the  wildlife  division. 
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